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Abstract 

For a set P of n points in M?, the Euclidean 2-center problem computes a pair of 
congruent disks of the minimal radius that cover P. We extend this to the (2, /c) -center 
problem where we compute the minimal radius pair of congruent disks to cover n — k 
points of P. We present a randomized algorithm with 0{nk^ log^ n) expected running 
time for the (2, /c)-center problem. We also study the (p, /c)-center problem in under 
the -metric for p = {4, 5}. We propose an k^'^^^nlog n algorithm for computing a ioo 
(4, /c) -center and an k^^^^n log^ n algorithm for computing a £oo (5, /c) -center. 

1 Introduction 

Let P be a set of n points in R^. For a pair of integers Q < k < n and p > 1, a family of p 
congruent disks is called a (p, k) -center if the disks cover at least n — k points of P; (p, 0)- 
center is the standard p-center. The Euclidean (p, A:) -center problems asks for computing a 
(p, A:)-center of P of the smallest radius. In this paper we study the (2, A:)-center problem. We 
also study the (p, A:)-center problem under the ^oo -metric for small values of p and k. Here we 
wish to cover all but k points of P by p congruent axis-aligned squares of the smallest side 
length. Our goal is to develop algorithms whose running time is n{k log n)^^^\ 

Related work. There has been extensive work on the p-center problem in algorithms and 
operations research communities fH [131 EOl [HI- If p is part of the input, the problem is NP- 
hard [24J even for the Euclidean case in M^. The Euclidean 1 -center problem is known to 
be LP-type |[22ll . and therefore can be solved in linear time for any fixed dimension. The 
Euclidean 2-center problem is not LP-type. Agarwal and Sharir ||2| proposed an 0{in? log^ n) 
time algorithm for the 2-center problem. The running time was improved to 0{n log^^^^ n) by 
Sharir [26|. The exponent of the logn factor was subsequently improved in ||T4l[5l. The best 
known deterministic algorithm takes 0{n log^ n log^ log n) time in the worst case, and the best 
known randomized algorithm takes 0{n log^ n) expected time. 

There is little work on the (p, A:) -center problem. Using a framework described by Ma- 
tousek [21 J, LP-type problems, with k violations and basis size 3, can be solved in 0{n log k + 
k^rf) time, for any e > 0. This is improved by Chan (61 to 0(n(3{n) logn + fc^n^) expected 
time, where /?(•) is a slow-growing inverse-Ackermann-like function and ^ > 0. The (1, k)- 
center problem is LP-type with basis size 3, so these bounds apply. Matousek |[2T1l also gives a 
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more general results for LP-type problems with k violations and with basis size c that runs in 
0{nk^) time, if it is a feasible case where a solution with no violations exists. In the infeasible 
case, no solution exists without violations and the algorithm runs in 0{nk'^^^) time. In fact, he 
shows in the feasible (resp. infeasible) case that there are 0{k^) (resp. 0{k^^^)) basis with at 
most k violations and his algorithm visits all of them by a path of length 0{k^) (resp. 0{k^^^)) 
where consecutive basis in the path differ by inserting or deleting one constraint. 

The p-center problem under ^oo -metric is dramatically simpler. Sharir and Welzl fT7\ show 
how to compute the ioo p-center in near- linear time for p < 5. In fact, they show that the 
rectilinear 2- and 3-center problems are LP-type problems and can be solved in 0{n) time. 
Also, they show the 1 -dimensional version of the problem is an LP-type problem for any p, 
with combinatorial dimension 0{p). Thus applying Matousek's framework 121 L the ioo (p, k)- 
center in M? for p < 3, can be found in 0(fc^(^)n) time and in 0{k^^^^n), for any if the 
points lie in R^. 

Our results. Our main result is a randomized algorithm for the Euclidean (2, fc)-center prob- 
lem in 'E? whose expected running time is 0{nk^ log^ n). We follow the general framework of 
Sharir and subsequent improvements by Eppstein. We first prove, in Section 2, a few structural 
properties of levels in an arrangement of unit disks, which are of independent interest. 

As in EE nil, our solution breaks the (2, A:)-center problem into two cases depending on 
the distance between the centers of the optimal disks; (i) the centers are further apart than the 
optimal radius, and (ii) they are closer than their radius. The first subproblem, which we refer 
to as the well-separated case and describe in Section 3, takes 0{k^n log^ n) time in the worst 
case and uses parametric search ll23l . The second subproblem, which we refer to as the nearly 
concentric case and describe in Section 4, takes 0(A:''nlog^ n) expected time. Thus we solve 
the (2, A:)-center problem in 0{k^n log^ n) expected time. We can solve the nearly concentric 
case and hence the (2, A:)-center problem in 0{k^n^^^) deterministic time, for any S > 0. We 
present near-linear algorithms for the ioo [p^ /c)-center in for p = 4, 5. The ^oo (4, /c)-center 
problem takes 0{k^^^^n\og n) time, and the ^oo (5, A:)-center problem takes 0{k^^^^n\o^ n) 
time. We have not made any attempt to minimize the exponent of k. We believe that it can be 
improved by a more careful analysis. 



2 Arrangement of Unit Disks 

Let D = {Di, . . . , Dn] be a set of n unit disks in R^. Let A{V) be the arrangement of 
A{V) consists of 0{n?) vertices, edges, and faces. For a subset !3? C 2), let = Cloe^R-^ 
denote the intersection of disks in Each disk in 3? contributes at most one edge in 
We refer to 3{X) as a unit-disk polygon and a connected portion of 83 {^) as a unit-disk curve. 
We introduce the notion of a level in yi(D), prove a few structural properties of levels, and 
describe a procedure that will be useful for our overall algorithm. 



^The arrangement of T> is the planar decomposition induced by T)\ its vertices are the intersection points of 
boundaries of two disks, its edges are the maximal portions of disk boundaries that do not contain a vertex, and its 
faces are the maximal connected regions of the plane that do not intersect the boundary of any disk. 
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Figure 1 : (a) yi(D) ; shaded region is yi<i (D) ; filled (resp. hollow) vertices are convex (resp. concave) 
vertices of yi<i(D); covering of yi<i(D) edges by six unit-disk curves, (b) yi(r+); shaded region is 
yi<i(r+); and the covering of yi<i(r+) edges by two concave chains, (c) yi(r~); shaded region is 
yi<i (r~); and the covering of yi<i (r~) edges by two convex chains. 



Levels and their structural properties. For a point x G M?, the level of x with respect to 
D, denoted by A(x, !D), is the number of disks in D that do not contain x. (Our definition of 
level is different from the more common definition in which it is defined as the number of disks 
whose interiors contain x.) All points lying on an edge or face of A{V) have the same level, 
which we denote by A(0). For k < n, let Ak{T)) (resp. yi</e(D)) denote the set of points in 
whose level is k (resp. at most A:); see Figure[l} By definition, Aq{'D) = A<q{T}) = 
The boundary of A<k{'^) is composed of the edges of A{V). Let v G dDi D dD2, for 
Di^D2 G D, be a vertex of 5yi</e(D). We call v convex (resp. concave) if yi</e(D) lies in 
DinD2 (resp. D1UD2) in a sufficiently small neighborhood of v; see Figure[T];a). dA<o{T)) is 
composed of convex vertices. We define the complexity of A^ki'D) to be the number of edges 
of A(V) whose levels are at most k. Since the complexity of A<o{T)) is n, the following 
lemma follows from the result by Clarkson and Shor |[ TO l (see also Sharir ||25]| and Chan LZI). 

Lemma 2.1. [TTOl For k > 0, the complexity ofA^ki"^) is 0{nk). 

Remark. The argument by Clarkson and Shor can also be used to prove that yi</e(D) has 
0{k^) connected components and that it has 0(A:^) local minima in (+2/)-direction. See also 
llQlEIl- These bounds are tight in the worst case; see Figure [2j 

It is well known that the edges in the <A;-level of a line arrangement can be covered by A: + 1 
concave chains ifTSll , as used in IIllO. We prove a similar result for yi<jt(D); it can be covered 
by 0{k) unit-disk curves. 

For a disk Di, let 7^^ (resp. 7^") denote the set of points that lie in or below (resp. above) Di\ 
d^^ consists of the upper semicircle of dDi plus two vertical downward rays emanating from 
the left and right endpoints of the semicircle — we refer to these rays as left and right rays. 
The curve d^^ has a similar structure. See Figures [ifb) and (c). Set r+ = {7+ | 1 < i < n} 
and r~ = {7^~ | 1 < ^ < n}. Assuming that the x-coordinates of the centers of all disks in 2) 
are distinct, each pair of curves d^^ ^ djj' intersect in at most one point. (If we assume that the 
left and right rays are not vertical but have very large positive and negative slopes, respectively, 
then each pair of boundary curves intersects in exactly one point.) We define the level of a 
point with respect to r+, F", or r+ U F" in the same way as with respect to D. A point lies in 
a disk Di if and only if it lies in both 7^^ and 7^", so we obtain the following inequalities: 



max{A(x, r+), A(x, F")} < A(x, D). 



(2.1) 
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A(x, V) < A(x, r+ U r-) < 2A(x, V). (2.2) 

We cover the edges of A<k{^^) by concave chains as follows. The level of the (A:+l)st 
rightmost left ray is at most /c at y = — oc. Let pi be such a ray, belonging to j^. We trace 
dj^, beginning from the point at = — oc on pi, as long as dj^ remains in yi<jt(r+). We 
stop when we have reached a vertex v G yi</e(r+) at which it leaves yi</e(r+); is a convex 
vertex on A<k{^^)' Suppose v = dj^ D djj~. Then dA<k{^^) follows ^7+ immediately 
to the right of v, so we switch to d^j~ and repeat the same process. It can be checked that we 
finally reach y = — oc on a right ray. Since we always switch the curve on a convex vertex, 
the chain we trace is a concave chain composed of a left ray, followed by a unit-disk curve 
and then followed by a right ray. Let Aq", A^, . . . , A^ be the /c + 1 chains traversed by 
this procedure. These chains cover all edges of yi</e(r+), and each edge lies exactly on one 
chain. Similarly we cover the edges of yi</e(r~) by A: + 1 convex curves Aq , A^, . . . , A^. 
Let S = {^q", . . . , ^q", . . . , be the family of unit-disk curves induced by these convex 



and concave chains. By (2.1 ), S covers all edges of yi<jt(D). Since a unit circle intersects a 



unit-disk curve in at most two points, we conclude the following. 

Lemma 2.2. The edges ofA^kC^) be covered by at most 2/c + 2 unit-disk curves, and a 
unit circle intersects 0(k) edges ofA^kC^)- 

The curves in S may contain edges of A{T)) whose levels are greater that k. If we wish to 
find a family of unit-disk curves whose union is the set of edges in yi<jt(!D), we proceed as 
follows. We add the x-extremal points of each disk as vertices of A{T)), so each edge is now 



x-monotone and lies in a lower or an upper semicircle. By (2.1 ), only 0{k) such vertices lie 
in yi</e(D). We call a vertex of yi</e(D) extremal if it is an x-extremal point on a disk or an 
intersection point of a lower and an upper semicircle. An extremal vertex of the latter type 
is an intersection point of G S. Since each such pair intersects in at most two points, 

there are 0(A:^) extremal vertices. For each extremal vertex v we do the following. If there is 
an edge e of yi<jt(D) lying to the right of v, we follow the arc containing e until we reach an 
extremal vertex or we leave yi</e(D). In the former case we stop. In the latter case we are at 
a convex vertex of 5yi</e(2)), and we switch to the other arc incident on and continue. 
These curves have been drawn in Figure [TJa). This procedure returns an x-monotone unit-disk 
curve that lies in yi<jt(!D). It can be shown that this procedure covers all edges of yi</e(D). If 
yi<jt(D) is represented as a planar graph, we can compute these curves in time proportional to 
the number of edges in A^ki"^)- We thus obtain the following: 

Lemma 2.3. Let T) be a set ofn unit disks in M?. Given yi</e(2)), we can compute, in time 
0{nk), a family of 0{k^) x-monotone unit-disk curves whose union is the set of edges of 

Remark. Since A^ki"^) can consist of connected components, the 0(/c^) bound is 

tight in the worst case; see Figure |2] 

Dynamic Data Structures. We need a dynamic data structure for storing a set D of unit 
disks that supports the following two operations: 

• (Ol) Insert a disk into D or delete a disk from D; 
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Figure 2: Lower bound. yi<2(2)) (shaded region) has 4 connected components. The right 
image is zoomed in of the center of the left image. 

• (02) For a given k, determine whether yi</e(2)) 7^ 0. 

Hershberger and Suri |[T9ll , describe how to maintain J(D) under insertion/deletion in 0(log n) 
time per update and how to find the point in U(!D) with the smallest ^/-coordinate in 0(log n) 
time. We use this in conjunction with Matousek's algorithm f2T| for visiting all basis of an 
LP-type problem with at most k violations. Specifically we examine the LP-type problem of 
finding the smallest ^/-coordinate of U(D) with k violations, which has a basis size of 2 and 
may be infeasible. Thus the path to visit all basis is of length 0{k^) and using Hershberger and 
Suri's data structure we traverse each step of the path in 0(log n) time by inserting or deleting 
a constraint and finding the discs defining the minimal ^-coordinate. 

Lemma 2.4. There exists a dynamic data structure for storing a set ofn unit disks so that (Ol) 
can be performed in 0(log n) time, and (02) takes 0{k^ log n) time. 

Agarwal and Matousek [11 provide a data structure that can maintain the value of the radius 
of the smallest enclosing disk under insertions and deletions in O(n^) time per update, for any 
5 > 0. We combine this with Matousek's algorithm for LP-type problems, specifically for the 
(1, A:) -center problem. Similar to the above data structure, the algorithm determines a path of 
length 0(/c^) to traverse all basis with at most k violations, and each is traversed in 0{n^) time 
by handling an insertion or deletion using Agarwal and Matousek's data structure. 

Lemma 2.5. The exists a dynamic data structure for a set ofn points such that under inser- 
tion/deletion of a pointy it can return the answer to the (1, k)-center problem in 0{k^n^), for 
any S > 0. 

3 Well-Separated Disks 

In this section we describe an algorithm for the case in which the two disks Di, D2 of the 
optimal solution are well separated. That is, let ci and 02 be the centers of Di and D2, and let 
r* be their radius. Then ||ciC2|| > r*; see Figure [s] Without loss of generality, let us assume 
that ci lies to the left of 02. Let be the semidisk lying to the left of the line passing through 
ci in direction normal to C1C2. A line i is called a separator line if i^i H = and i separates 
from D2, or Di D D2 9 and i separates from the intersection points dDi H dD2. 
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Figure 3: Let ^ is a separator line for disks Di and D2. 



We first show that we can quickly compute a set of 0(A:^) lines that contains a separator line. 
Next, we describe a decision algorithm, and then we describe the algorithm for computing Di 
and D2 provided they are well separated. 



Computing separator lines. We fix a sufficiently large constant h and choose a set ?7 = 
{ui^ . . . , C of directions, where Ui = (cos(27rz//i), sin(27rz//i)). 

For a point p G and a direction Ui, let be the projection of p in the direction normal 
to Ui. Let pW = (p]^' , . . . , ) be the sorted sequence of projections of points in the direction 
normal to Ui. For each pair a, b such that a + 6 < A:, we choose the interval 6^^\ = [p^a , P^n-b\ 



and we place 0(1) equidistant points in this interval. See Figure 



3;a). LetL^^bethe set of 



(oriented) lines in the direction normal to Ui and passing though these points. Set 

i= U 4V 

l<i<h 
a+b<k 

The set L can be computed in 0(A:^nlogn) time. We claim that L contains at least one 
separator line. Let Ui e U ho the direction closest to C1C2. Suppose pa and pn-b are the first 
and the last points of P in the direction Ui that lie inside D1UD2. Since \P\ (Di U 1^2) | < k, 
a-\-h <k. Let be the extreme point of I? ^ in direction ui and let q2 be the extreme point of 
D2 \ Di in direction —ui. Since ui is within a small constant angle of C1C2 

(^2 - qi.Ui) > a{q2 - qi,CiC2) = -(C2 - Ci,CiC2) > -{Pn-b -Pa.Ui), 

z o 

where a < 1 is a constant depending on h. Hence if at least Q/a points are chosen in the 
interval 5^^^, then one of the lines in l}-^^ is a separator line. We conclude the following. 

Lemma 3.1. We can compute in 0{k^nlogn) time a set L of 0{k^) lines that contains a 
separator line. 
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Let Di,D2 be a (2, A:) -center of P, let ^ G L be a line, and let P~ C P be the set of points 
that lie in the left halfplane bounded by i. We call a (2, fc) -center consistent with i if 

P~ n {Di U D2) C Di, the center of Di lies to the left of i, and dDi contains at least one 
point of P~. We first describe a decision algorithm that determines whether there is a (2, k)- 
center of unit radius that is consistent with i. Next, we describe an algorithm for computing 
a (2, A:)-center consistent with i, which will lead to computing an optimal (2, A:)-center of P, 
provided there is a well- separated optimal (2, fc)-center of P. 



Decision algorithm. Let ^ G L be a line. We describe an algorithm for determining whether 
there is a unit radius (2, A:)-center of P that is consistent with i. Let P~ (resp. P+) be the 
subset of points in P that lie in the left (resp. right) halfplane bounded by i; set n~ = |. 
n+ = l^^l- Suppose P^i, is a unit-radius (2, A:)-center of P consistent with i, and let ci, C2 
be their centers. Then P" n (Di U D2) C P)i and \P~ n Di\ >n~-k. For a subset Q C P, 
let X)((5) = \ q ^Q} where P)(g) is the unit disk centered at q. Let D" = and 

X)+ = D(P+). For a point x G let D+ = {P> G D+ | x G P>}. Since aP>i contains a 
point of P~ and at most k points of P~ do not lie in Di, c\ lies on an edge of A<^]^{^~). 
We first compute yi</e(D~) in 0(n/clogn) time.For each disk D G D^, we compute the 



intersection points of dD with the edges of yi</e(D ). By Lemma 2.2 there are 0{nk) such 



intersection points, and these intersection points split each edge into edgelets. The total number 



of edgelets is also 0{nk). Using Lemma 2.2 we can compute all edgelets in time 0{nk log n), 
because each disk boundary from !D+ intersects at most 0{k) edges of A^k{'^~) and each 
intersection can be found in O(logn) time be examining the covering unit disk curves. All 
points on an edgelet 7 lie in the same subset of disks of D+, which we denote by 2)+. Let 
P+ C P+ be the set of centers of disks in D+, and let = A(7, D~) be the level of 7 
in D~. A unit disk centered at a point on 7 contains P+ and all but points of P~. If at 
least k' ^ k — points of P+ \ P+ can be covered by a unit disk, which is equivalent to 
A<k' (2^^ \ 2^7) being nonempty, then all but k points of P can be covered by two unit disks. 

When we move from one edgelet 7 of yi</e(!D~) to an adjacent one 7' with a as their 
common endpoint, then X)+ = 2)+ (if a is a vertex of yi<A;(X>")), = X>| U {D} (if 
a e dD and 7' C {D}), or D+ = D| \ (if a G aP> and 7 C D). We therefore 
traverse the graph induced by the edgelets of yi</e(D) and maintain in the dynamic data 
structure described in Section [2] as we visit the edgelets 7 of yi</e(D~). At each step we 
process an edgelet 7, insert or delete a disk into 2)+, and test whether A<j{V^) = where 
j = k — A(7, D~). If the answer is yes at any step, we stop. We spend 0{k^ logn) time at 
each step, by Lemma [2!4} Since the number of edgelets is 0{nk), we obtain the following. 



Lemma 3.2. Let P be a set of n points in M?, £ a line in L, and < k < n an integer. We 
can determine in 0{nk'^ logn) time whether there is a unit-radius (2, k) -center of P that is 
consistent with i. 



Optimization algorithm. Let ^ be a line in L. Let r* be the smallest radius of a (2, A;)- 
center of P that is consistent with i. Our goal is to compute a (2, A:)-center of P of radius 
r* that is consistent with i. We use the parametric search technique Ii23ll — we simulate the 
decision algorithm generically at r* and use the decision algorithm to resolve each comparison, 
which will be of the form: given tq G is ro < r*? We simulate a parallel version of the 
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decision procedure to reduce the number of times the decision algorithm is invoked. Note that 
we need to parallelize only those steps of the simulation that depend on r*, i.e., that require 
comparing a value with r*. Instead of simulating the entire decision algorithm, as in [14J, we 
stop the simulation after computing the edgelets and return the smallest (2, A:)-center found so 
far, i.e., the smallest radius for which the decision algorithm returned "yes." Since we stop the 
simulation earlier, we do not guarantee that we find the a (2, A:)-center of P of radius r* that is 
consistent with i. However, as argued below this is sufficient for our purpose. 

Let P~, P+ be the same as in the decision algorithm. Let D~, etc. be the same as above 
except that each disk is of radius r* (recall that we do not know the value of r*). We simulate 
the algorithm to compute the edgelets of A<k{'D~) as follows. First, we compute the <k^^ 
order farthest point Voronoi diagram of P~ in time 0{n log n + nk'^) [4J. Let e be an edge of 
the diagram with points p and q of P~ as its neighbors, i.e., e is a portion of the bisector of p 
and q. Then for each point x G e, the disk of radius | \xp\ \ centered at x contains at least n~ — k 
points of P~. We associate an interval Se = | x G e}. By definition, e corresponds 

to a vertex of yi<jt(D~) if and only if r* G Se', namely, if \\xp\\ = r*, for some x G e, then 
X is a vertex of yi</e(2)~), incident upon the edges that are portions of dD{p) and dD{q). 
Let X be the sorted sequence of the endpoints of the intervals. By doing a binary search on 
X and using the decision procedure at each step, we can find two consecutive endpoints in 
X between which r* lies. We can now compute all edges e of the Voronoi diagram such that 
r* G 5e. We thus compute all vertices of yi<jt(!D~). Since we do not know r*, we do not have 
actual coordinates of the vertices. We represent each vertex as a pair of points. Similarly, each 
edge is represented as a point p G P~, indiciating that e lies in dD{p), and it can be computed 
using the cells of the Voronoi diagram. Given a vertex of yi</e(D~) and an outgoing edge, 
represented by the point p G P~, we can compute the other endpoint as the next edge of 
the Voronoi cell of the p that is a point in A<{V~) by walking around the boundary of the 
cell. Once we have all the edges of yi</e(P~), we can construct the graph induced by them 
and compute 0(A:^) x-monotone unit-disk curves whose union is the set of edges in yi</e(P~), 
using Lemma [23] Since this step does not depend on the value of r*, we need not parallelize 
it. Let S = {^^, . . . , ^u}, u = 0(k'^), be the set of these curves. 

Next, for each disk D G and for each G S, we compute the edges of that dD 
intersects, using a binary search. We perform these 0{nk^) binary searches in parallel and use 
the decision algorithm at each step. Incorporating Cole's technique lITTIl in the binary search, 
the decision procedure is invoked only O(logn) times. For an edge e G yi</e(D), let G T) 
be the set of disks whose boundaries intersect e. We sort the disks in !D+ by the order in which 
their boundaries intersect e. By doing this in parallel for all edges and using a parallel sorting 
algorithm for each edge, we can perform this step by invoking the decision algorithm 0(log n) 
times. The total time spent is 0{nk^ log^ n). 



Putting pieces together. We repeat the optimization algorithm for all lines in L and return 
the smallest (2, A:) -center that is consistent with a line in L. Since Lemma 3.1 shows that as 
long as the solution is well separated at least one line in L is a separator line for the optimal 
(2, A:)-center of P, the smallest radius returned must be that of the optimal (2, A:)-center of P. 
Hence, we conclude the following: 

Lemma 3.3. Let P be a set of n points in and < k < n an integer. If an optimal 
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(2, k)-center of P is well separated, then the (2, k)-center problem for P can be solved in 
0{nk^ log^ n) time. 



4 Nearly Concentric Disks 




Figure 4: Two unit disks Di and D2 or radius r* with centers closer than a distance r*. 

In this section we describe an algorithm for the case in which the two disks Di and D2 of 
the optimal solution are not well separated. More specifically, let ci and C2 be the centers of 
Di and D2 and let r* be their radius. This section handles the case where | |ciC2 1 1 < r*. 

First, we find an intersector point z of Di and D2 — a point that lies in DiH D2. We show 
how z defines a set 7 of 0{in?) possible partitions of P into two subsets, such that for one 
partition P^j, P \ Pij the following holds: {Di U L>2) H P = (P>i H Pij) U {D2 H (P \ Pij)). 
Finally, we show how to search through the set 7 in 0{k^n^^^) time, deterministically, for any 
5 > 0, or in 0{k^n log^ n) expected time. 

Finding an intersector point. Let C be the circumcircle of P n {Di U D2). Eppstein |[T4ll 
shows that we can select 0(1) points inside C such that at least one, z, lies in Di H D2. We 
can hence prove the following. 

Lemma 4.1. Let P be a set ofn points in M?. We can generate in 0{nk^) time a set Z of 
0{k^) points such that for any nearly concentric (2, k)-center Di, D2, one of the points in Z 
is their intersector point. 

Proof. Using Matousek's (2V\ algorithm for solving LP-type problems with violations, in 
0{k^n) time we can find the smallest circle that contains n — k points of P. Briefly, the 
algorithm runs by finding the three points defining the circumcircle, removing each one in 
turn, and recursing until k points have been removed. Matousek shows that if we keep track of 
which nodes in the recursion we reach and halt the recursion if we have seen that node before, 
then the size of the recursion tree is only 0(A:^). In the running of this algorithm we generate 
all circles which include exactly n — j points of P for < j < k. We claim that one of these 
circles must be C. 
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If the initial circle is not C, then it must have at least one point on its boundary which is not 
in P n {Di U D2). At least one path of the recursion removes this point. Since we can reach 
the point set P H {Di U £^2) in at most k steps, some step in this recursion must return C. 

Finally, since the area of Di U i?2 is a constant fraction of C when I?i,i?2 are nearly 
concentric, then by selecting a constant number of points in C one can be guaranteed to be an 
intersector point. □ 

Let z be an intersector point of Di and D2, and let p~ be the two rays from z to the 
points of dDindD2. Since Di and D2 are nearly concentric, the angle between them is at least 
some constant 9. We choose a set (7 C 5^ of /i = l^^^/O] uniformly distributed directions. 
For at least one u e U, the line i in direction u and passing through z separates p+ and p~, see 
Figure |4j We fix a pair z,u in Z x U and compute a (2, A:)-center of P, as described below. 
We repeat this algorithm for every pair. If Di and D2 are nearly concentric, then our algorithm 
returns an optimal (2, fc) -center. 

Fixing z and u. For a subset X C P and for an integer t > 0, let r^{X) denote the minimum 
radius of a (1, t)-center of X. Let P+ (resp. P~) be the subset of P lying above (resp. below) 
the X-axis; set n+ = |P+| and n~ = |^~|- Sort P+ = (p^, . . . in clockwise order 

and P~ = (p^, . . . in counterclockwise order. For < i < n+, < j < n~, let 

PiJ = {Pi^ • • • ^vt^Pi^ ' ' • ^Pj} and Qij = P \ Pij. For < t < fc, let 

m^,^. = max{r^(P,,,),r^-^(Q.,,)}. 

For < t < fc, we define an n+ x n~ matrix such that M^(z, j) = m^^. 

Suppose z is an intersector point of Di and P)2, ^ separates p+ and and p+ (resp. p~) 
lies between (resp. p^^p^^^l Then P n (P>i U P>2) = {Pa^b H P>i) U (Q^,^ U P>2); 
see Fig|4j If \Pa^b \ ^ t, then r* = ^. The problem thus reduces to computing 

/x(z, u) — minm^ • 

where the minimum is taken over < z < n+, < j < n~, and < t < k. For each t, we 
compute iJ{z, u) = min^^j m\ - and choose the smallest among them. 

We note two properties of the matrix that will help search for //^(z, u): 

• (PI) If r^{Pij) > r^~'^{Qij) then mlj < m^, j; for i' > i and > j. These partitions 
only add points to Pij and removes points from Qij, and thus cannot decrease r^{Pij) 
or increase r^~'^{Qij). Similarly, if r^~'^{Qij) > r^{Pij), then mlj < m\, ■, for i' <i 
and/<j. 

• (P2) Given a value r, if r^{Pij) > r, then m^,^, > r for > z and > j, and if 
r^{Qij) > r, then m^, j, > rfovi' <i and < j. 

Deterministic solution . We now have the machinery to use a technique of Frederickson and 
Johnson ifTSll . For simplicity, let us assume that n+ = n~ = 2^^^ where r = [log2 n \ + 0(1). 
The algorithm works in r phases. In the beginning of the /ith phase we have a collection JAh 
of 0(2^) submatrices of M^ each of size (2^-^+^ + 1) x (2^-^+^ + 1). Initially Mi = 
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{M*}. In the hth phase we divide each matrix A'' e Mh into four submatrices each of size 
(2^-'* + 1) X (2^-'* + 1) that overlap along one row and one column. We call the cell common 
to all four submatrices the center cell of N. Let JA'^ be the resulting set of matrices. Let 
C = {(n, ji ),•••, (^s, js)} be the set of center cells of matrices in M/^. We compute Tn\^ -^ for 
each 1 < / < 5. We use (PI) to remove the matrices of M/^ that are guaranteed not to contain 
the value //^(z, u). In particular, if = r^{P^^j^) and there is a matrix N E JA'^ with the 
upper-left corner cell (i^ ^^^h that i' < ii and f < ji, then we can remove N. Similarly if 
31 ^ ^^~^(Qhj) there is a matrix N G M'^ with the lower-right corner cell j') such 
that i' > ii and j' > ji, we can delete N. We then set M'^ to M/^+i. 




Figure 5: Example of running deterministic algorithm through 3 phases. Shaded regions have 
been pruned. Center cells are darkened. 



Lemma 4.2. Before the hth phase consider a diagonal from large i and j to small i and j that 
passes through at least one center cell of a matrix N G M/^. It passes through at most one 
more center cell of a matrix N' G M/^. 
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Proof. We show this inductively. The base case is clearly true for the single center cell in Mi. 
Assume it is true for M/^, then we show it is true for M/^+i. See Figure |5j We consider two 
cases, first the diagonal passes through a center cell of M/^. In this case if it passes through two 
center cells of M/^, then it passes through 4 center cells of M'^, but the pruning step eliminates 
at least two of them. In the second case, the diagonal does not pass through a center cell of 
M/j. We can bound the number of center cells of matrices it passes through in JA'^ to 4 using 
the inductive hypothesis. Consider one of the interior center cells {i^j) ^ N ^ M'^ it passes 
through, neither the first not the last. When the pruning step for the matrix in Mh that contains 
N is called, it either eliminated the other matrixes in M'^ that the diagonal passes to before 
or after N. If the diagonal passes through 3 center cells in M^^, then this reduces it to two, if 
the diagonal passes through 4 center cells, the applying this analysis to both interior matrices 
reduces it to two. □ 

Lemma [4~2| implie s that 0{n) cells remain in M'^ after the pruning step and that they can be 
connected by two monotone paths in M^ which consists of 0(n) cells. Since Pij differs from 
Pi-ij and Pij-i by one point, we can compute rnj^ j^ for all (i/, j/) G C using Lemma 
total time 0{k^n^^^). Hence, each phase of the algorithm takes 0{k^n^^^) time. 



2.5 



m 



Lemma 4.3. Given z ^ Z, u ^ U, and < t < k, iJ{z^u) can be computed in time 
0{k^n^+^),forany6> 0. 



Randomized solution. We can slightly improve the dependence on n by using the dynamic 
data structure in Section [2] and (P2). As before, in the hih phase, for some constant c > 1, we 
maintain a set Mh of at most c2^ submatrices of each of side length 2^^"^+^ + 1, and their 
center cells C. Each submatrix is divided into four submatrices of side length 2^~^ + 1, forming 
a set M'^. To prune M^^, we choose a random center cell (z, j) from C and evaluate r = ml j in 
0{k^n) time. For each other center cell (z^, f) G C, m\, ■, > r with probability 1/2, and using 
(P2), we can remove a submatrix from M^^. More specifically, if m^, j, > r, then any matrix 
N G M'^ with an lower right corner (z^, j') such that i' < i and j' < j or a upper left corner 
{i" ^j") such that i" > i and j^^ > j, then we can prune N from M^^. Eppstein |[T4l proves that 
by repeating this process a constant number of times, we expect to reduce the size of M'^ to 
c2^+i. 

On each iteration we use the dynamic data structure described in Section |2j For 0(n) inser- 
tions and deletions, it can compare each center cell from C to r in 0{k^n log^ n) time. Thus, 
finding u) takes expected 0{nk^ log^ n) time. 

Lemma 4.4. Given z ^ Z, u ^ U, and < t < k, fi^(z^u) can be computed in expected time 
0{k^ log^ n). 



Putting pieces together. By repeating either above algorithm for all < t < A: and for 

all pair (z, u) ^ Z xU.wq can compute a (2, A:)-center of P that is optimal if Di and D2 are 
nearly concentric. Combining this with Lemma [33} we obtain the main result of the paper. 

Theorem 4.1. Given a set P ofn points in M? and an integer k > 0, an optimal (2, k)-center 
of P can be computed in 0{k^in}^^) (deterministic) time, for any 5 > {) or in 0{k^nlog^ n) 
expected time. 
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5 The {p,k)-Center Problem Under the i^o Metric 

This section focuses on the ioo version of the (p, A:)-center problem, and hence all references 
to the (p, fc) -center problem herein are referring to the ioo variant. We use extensively that for 
p < 3, the (p, 0) -center problem is LP-type f27|, and thus for p < 3 the (p, A:) -center problem 
can be solved in 0{k^^^^n) time. We also use that if all points lie in R^, then the (p, 0)-center 
problem is LP-type for any p > 0, with combinatorial dimension 0{p), and thus in R^, the 
(j), A:)-center problem can be solved in 0{k^^'^^n) time. 

Like in the ^2 variant, we first study the decision version of the dual problem; here an 
arrangement of unit squares. Let S = {S^i, . . . , 5^^} be a set of n unit squares (side length 
1) in R^. Let A{§) be the arrangement of S. We say a point q stabs a square 5 G S if g G 5. 
Let S(g) C S be the set of squares stabbed by q. 

We seek to determine whether there exists a placement of p points Q = {^1, • • • , ^p} such 
that UgGQ^(^) ^ ^ — k. We refer to this as the (p, fc) -stabbing decision problem. All of 
our algorithms also return a solution if one exists. By replacing each point in the (p, A:) -center 
problem with a unit square centered at that point, then the p stabbing points of the (p, k)- 
stabbing decision problem serve as the center points of unit squares that contain n — A; of the 
original point set. 



Structure. We start by reviewing structure observed by Sharir and Welzl ll27ll about the 
(p, 0) -stabbing decision problem. 

If a horizontal or vertical line i passes through all 5 G S, then this (p, 0)-center decision 
problem reduces to a variant in R^ because any stabbing point q can be replaced with q\ the 
closest point on i to q, so that §{q) < §{q'). We can then solve the (p, 0) -stabbing decision 
problem in 0{n) time or the (p, A:)-center problem in 0{nk^^P'^) time. We henceforth assume 
that this is not the case. 

Let describe the line passing through the right boundary of the leftmost square. Similarly, 
let (resp. i^, i^) describe the line passing through the left (resp. bottom, top) boundary of 
the rightmost (resp. topmost, bottommost) square. Let Hq describe the rectangle bounded on 
its left side by i^, its right side by i^, its bottom side by i^, and its top side by i^. (See Figure 
[6j) Ho must have positive area otherwise a horizontal or vertical line would pass through the 
set of all squares. 

Let Hq n , for X G {L, i?, T, B}, describe the four boundary segments of Hq. Call the 
intersection of two boundary segments a corner of Hq. If the (p, 0) -stabbing decision problem 
has a solution, we claim that each boundary segment of Hq contains a stabbing point in a 
solution of the (p, 0) -stabbing decision problem (in particular, the solution of p stabbing points 
contained in the smallest rectangle). For instance, if Hq H does not contain a stabbing point, 
then we can replace the point q stabbing the leftmost square with another point q^ on Hq D 
such that S(g) < §{q'). 

If a stabbing point q lies on corner, it lies on two boundary segments at once, and we can set 
= §\§{q) and then solve the — 1, A:) -stabbing decision problem on S^ Of course, we 
don't know which corner is a stabbing point, but there are a constant number and we can try 
them all. 

Define i\ (resp. I^t^^^b) as the line through the right (resp. left, bottom, top) boundary 
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of the j^th leftmost (resp. j^th rightmost, j^th topmost, j^th bottommost) square. We can 
also define the rectangle HjL jR jT jB which is defined by the intersection of halfspaces defined 
by lines ij^ ^ ^fn ^ ^Jt ^ and i^s • We actually want to be slightly careful since one square may be 
in the j^th leftmost and j^th topmost squares. We count squares first from left and right, then 
those remaining from top and bottom. Let SjL jR jT jB be the set of squares which intersect 

Dynamic data structure. We will need a data structure to be able to maintain Hq and 8 
under the removal of the set §{q) for a possible stabbing point q. Sharir and Welzl [27] provide 
a data structure that stores a set of canonical subsets, such that under this operation §\§{q) 
can be stored as the union of O(logn) (not necessarily disjoint) canonical subsets. The new 
boundary lines of Hq can be constructed in 0(log n) time from the 0(log n) subsets. 

The structure is built, and extended to handle outliers, as follows. In the x- and ^/-directions 
store binary trees of S sorted by their coordinates. Each node in the tree stores a canonical 
subset of all squares in its subtree. For a query point q, we can return all squares that cannot 
intersection q based on x- and ^/-coordinates independently, as a set of O(logn) canonical 
subsets each. The union is S \ S(g). We may need to build this data structure p — 1 levels 
deep on each canonical subset for solving the (p, A:) -stabbing decision problem. To construct 
Hq quickly, we can find the maximum and minimum square in x and y coordinate over all 
0(log n) canonical subsets. To instead construct HjL jRjT jB, we can find the minimum x 
coordinate in 0{j^ log n) time and similarly for j^, j^, and j^; thus constructing HjL jR jT jB 
can be done in 0{k log n) time, where j^^j^, j"^, < k. 

5.1 The (4, /c)-Stabbing Decision Problem 

Firstwechoosepositiveintegral values and such that j^+j^+j^+j^ < k-\-4 

and create HjL jR jT jB. If j^, j^, j^, and are chosen correctly, then HjL jR jT jB is the 
smallest rectangle that contains the 4 stabbing points. If the decision is true, then one of this set 
of 0(A:^) rectangles must match the solution because it can not exclude more than k rectangles 
in any one direction. In what follows, we assume we have chosen ^j^ correctly, but 

in the full algorithm we try each until we find a solution. If HjL jR jT jB has non positive area 
then we can solve the problem in R^. We then see if one of the corners, q, of HjL jR jT jB 
can be a stabbing point by solving the (3, k — {j^ + + + — 4))-stabbing decision 
problem on SjL jR jT jB \ §{q). If the answer is negative for each corner, and we assume that 
we have chosen j^, j^, j^, correctly, then each boundary segment of HjL jR jT jB must 
contain a distinct stabbing point. Let C SjL jR jT jB be the subset so that each 5 G 8^ does 
not intersection dHjL jR jT jB — these squares must be totally contained in HjL jR jT jB . Let 

In the following we assume that HjL jR jT jB is the smallest rectangle to contain all stabbing 
points and to simplify notation we set = SjL jR jT jB \ S^, = A: — (j^ + + + — 
4) — k^, and H = HjL jR jT jB. Finally, we assume that the solution to the (p, /^)-stabbing 
decision problem on has no point on the corners of H. 
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Figure 6: Structure of a 4-center problem with ^oo-distance. Rectangle H = i?i,2,3,2 is shaded 
and bounded by lines i^^i^^ ^i^ on the top, right, bottom, and left sides, respectively. The 
four centers appear on the four sides of H labeled. There are five outliers squares shown in 
bold. 



4 Rotating Calipers. We can now apply a rotating calipers type technique with four calipers, 
with one point on each edge of H. Since each square can intersect each edge of H at most 
twice, the boundary of H is divided into 0(n) regions such that all points within a region of 
the boundary intersect the same set of squares. Squares can intersect more than one edge of 

H, either by also containing a corner point (i.e. left and top), intersect two opposite sides (i.e. 
left and right), or both (i.e. left, top, and right). In the third case when a square intersects three 
sides it must entirely contain one of those sides, and thus any point chosen on that side must 
stab that square and we can ignore it. Also only one pair, w.l.o.g. top and bottom, can have 
squares intersecting both, otherwise both pairs of opposite sides are shorter than a distance 

I, and any square intersecting a pair of opposite sides must entirely contain one of the other 
sides. Assuming the top and bottom edges are longer than 1 (so no square can intersect both 
the left and right edge) we consider two cases: where the point on the top side is right of the 
point on the bottom side, and vice versa. We focus on the first case and handle the other one 
symmetrically. 

We treat the subset of squares 82 ^ which intersect the top and bottom edges separately 
from the subset = \ 82 of the ones that only intersection only one edge or two adjacent 
edges. Each square S ^ §e describes one interval on the curve defined by dH. Thus, given a 
placement of four stabbing points, one on each boundary side, the squares from S e which are 
not stabbed lie in one of four intervals of dH bounded by the stabbing points. In the optimal 
solution let there be unstabbed squares in between the and q^, squares between 
q^ and q^, squares between q^ and q^, and squares between q^ and q^. For any values 
z^, z^, z^, and we can determine if there is placement of the stabbing points on dH that 
has exactly those many unstabbed squares in the associated intervals. Given a placement of q^ 
in bottommost region of the right boundary edge, we can try to place q^ skipping squares, 
then place q^ skipping squares, and finally q^ skipping squares. If there are squares 
remaining it is successful. If it is not successful at any placement step, then we shift q^ to the 
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next region up on the right boundary and try shifting the other stabbing points to the next region 
in a counter-clockwise direction to satisfy the constraints. If all attempts are unsuccessful for 
all placements of on the right edge, then this choice of i^, and is incorrect. Since 

there are only 0(n) regions, and each stabbing point is in each region at most once, since they 
only move counter-clockwise, this takes 0{n) time. 

Once a solution for has been found, we attempt find a solution for 82- These squares 
can be sorted left to right and a successful stabbing will have unstabbed squares from 82 
left of q^, squares between and q^, and squares right of q^, for some nonnegative 
integers i^, i^, i^. We can now adjust q^ and q^ such that the sets §e{(1^) and §e{(1^) do not 
change. The boundaries of the squares from 82 divide the regions into intervals so that within a 
interval §2{q^) and §2(0^) do not change. After preprocessing to find the left boundary of the 
rightmost square in 82 and the right boundary of the leftmost square in 82, in 0{i^ + i^) time 
we check if we can place q^ and q^ to satisfy and i^. Quickly checking the constraint 
requires preprocessing on the intervals created by the sorted ordering of 82 so each region 
contains the number of points stabbed and the number of unstabbed squares to the right. Thus 
if q^ is in a region so that it stabs s squares and there are r squares to the right of the region 
that q^ is in, then there are r — 5 — unstabbed squares from 82 between q^ and q^. If 
r — s — i^ then we return true, if not we go back to dealing with 8 e and shift the stabbing 
points in counter-clockwise order. 

Although, we do not know the values of i^, i^, , i^, i^, i^, and we do know that 
+ + + + + + = thus there are only 0{k^) possible values. For each 
set of values, we require 0{n) time to handle 8^ and for each step 0{k) time to handle 82, 
after preprocessing. Let Too(n,^, k) be the required time for the algorithms described above 
to solve the (p, fc)-stabbing decision problem on n unit squares. 

Lemma 5.1. Too{n,A,k) = 0{k^{Too{n,?>,k) + k^^n) + n\ogn) orjustO{k^{Too{n,?>,k) + 
k^^n)) if the squares are presorted along the x- and y-axis. 

Theorem 5.L Too{n,4,k) = 0{k^^^^n + nlog n) or just Oik^^^'^n) if the squares are pre- 
sorted along the x- and y-axis. 

5.2 The (5, fc)-Stabbing Decision Problem 

We first construct 0{k^) rectangles H = dHjL jR jT jB as above. To simplify notation, also 
assume that k squares intersect or lie inside of H and that at least one center must lie on each 
side of H. We now have to consider 3 cases. 

First, one of the centers lies on a corner of H. In this case, we can try all corners, remove the 
squares that intersect that comer and apply the algorithm for p = 4 on the remaining squares. 

Second, all 5 of the centers lie on the rectangle H (not its interior), but none lie on a corner. 
Third, 4 centers lie on the boundary of H, but none lie on a corner and the fifth center lies in 
the interior of H. These cases are more complicated and requires the dynamic data structure 
described above. 

In the second and third case we choose non-negative integers ii through ig such that ii +^2 + 
^3 + ^4 + ^5 + ^6 + ^7 + ^8 + ^9 = k. These determine which points are outliers and not contained 
in the 5 centers. We choose 0{k^) sets of integers and complete the following for each set. 
For what follows we assume we have chosen the correct set. Each side of H has at least one 
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center, and one side has two. We perform the following, assuming each side, in turn, has two 
centers; w.l.o.g. let it be the right side. Guess that the left most interval on the bottom edge of 
H contains the center point, p^. Let S^s = S \ S(p^) be the set of squares that do not contain 
p^. Using the above 4-level dynamic data structure, obtain S^s and construct H' = 
on §pB. The bottommost ii remaining squares have been designated as outliers, not to contain 
any center point. Now either the bottom left or the bottom right corner of must contain a 
center point. Check each case by the following; w.l.o.g. assume its the bottom right corner, 
p^. Create S^s \ Sps(p^) using the dynamic data structure, and recalculate 1+^2, 1,1+^3- 
Again the 12 rightmost squares and is bottommost squares are designated outliers. Now again 
either the bottom right or bottom left corner of H must be a center point. Check either, remove 
Z4 and Z5 outliers, and proceed as before removing the squares contained in the third center. 
This process repeats once more, removing and outliers, and squares containing the fourth 
center. There is now one center left to place. We remove is and zg outliers and can easily check 
if the last center can contains all remaining squares. If it cannot then we update p^ by sliding it 
to the next interval on the bottom edge of H. We update our 4-level data structure in 0(log^ n) 
time. This repeats until either all squares can be stabbed by the last center, meaning the result 
is true, or all intervals on the bottom edge of H have been tried, meaning the result is false. 

Accounting for the 0{k^) possible outliers to create the initial rectangle H, and the 0{k^) 
sets of integers ii . . . zg, the final running time is 0{nk^'^ log^ n). 

The third case is very similar to the second case. We consider a case where p^, the center on 
the interior of H, is above either the point on the left side or the right side of H. If this is not 
true, we would perform the process symmetrically by guessing a center on the top side instead 
of the bottom side. We can remove squares containing the first two center points the same way 
as in the second case. When there are three center points remaining, we can still claim that one 
lies on the corner of H, but its not necessarily a bottom corner. This just requires a few more 
cases to check. It follows that this third case also takes 0{nk^'^ log^ n) time. 

5.3 The (4, k)- and (5, /c)-Center Problem 

We can solve the original primal problem of determining whether a set of p squares can contain 
all but k points from an n point set. To find the minimum side length of the squares for this to 
be true we can use a matrix searching technique of Frederickson and Johnson [[T5l[T6l[T7l with 
0(log n) iterations of the above algorithm. The minimal side length of a square is necessarily 
the difference in x-coordinates between two points or the difference in ^-coordinates between 
two points. We implicitly store these two orderings along the columns of two matrices, X and 
y , corresponding to the x- and the ^-coordinates of the points. The cells contain the differences 
in their values, but are only computed as needed. Using monotone properties of these matrices 
we can search for the minimum such difference where our algorithm returns true. We take the 
minimum from both matrices. 

Theorem 5.2. Given a set Pofn points in and an integer k > 0, an optimal (4, k) -center 
ofP can be computed under the £oo-metric in ©(/c^^^^n log n) time. 

Theorem 5.3. Given a set P ofn points in and an integer k > 0, an optimal (5, k)-center 
ofP can be computed under the £oo-metric in 0{k^^^'^nlog^ n) time. 
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